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1.  INTRODUCTION 

In  the  past,  dynamic  behavior  of  periodically  laminated  medium  has  been  studied 
extensively.  A  review  of  the  literature  on  ex2u:t  zind  approximate  analyses  of  this 
problem  can  be  found  in  [1,2]. 

The  corresponding  problem  of  a  plate  of  finite  thickness  having  a  large  number 
of  periodic  anisotropic  layers  has  not  received  much  attention.  In  [2]  a  stiRness 
method  was  presented  for  studying  harmonic  wave  propagation  in  a  periodically  lam¬ 
inated  infinite  medixim.  This  method  is  also  well  suited  for  analyzing  a  finite  thick¬ 
ness  plate  with  many  layers.  In  this  approach  each  lamina  is  divided  into  several 
sublayers  and  the  displacement  distribution  through  the  thickness  of  each  sublayer  is 
approximated  by  polynomial  interpolation  functions  in  such  a  way  that  displace¬ 
ments  and  tractions  are  continuous  across  the  interfaces  between  adjacent  sublayers. 
Details  of  the  method  can  be  found  in  [2,3|.  Here  we  summarize  the  pertinent  equa¬ 
tions  and  present  numerical  results  showing  the  effect  of  soft  interface  layers  between 
reinforced  laminae  on  the  dispersion  of  guided  waves. 

The  particular  systems  considered  are  a  single  ply  (0  * )  and  a  cross-ply 
(0  *  /90  *  /^  * )  laminated  plate.  Each  lamina  is  assumed  to  be  reinforced  by  continu¬ 
ous  aligned  fibers,  so  that  it  can  be  modeled  as  a  transversely  botropic  medium  with 
the  symmetry  axis  aligned  with  the  fiber  direction.  The  geometry  of  a  typical  cross- 
ply  laminate  is  shown  in  Fig.  1,  which  also  shows  thin  interface  layers  between  adja¬ 
cent  laminae.  We  have  also  considered  the  case  when  all  the  fibers  are  in  the  same 
direction  (not  shown).  Although  it  is  possible  to  derive  an  exact  dispersion  equation 
for  propagation  of  harmonic  waves  in  the  plane  of  the  plate,  the  equation  is  rather 
complicated  when  the  number  of  layers  is  large.  This  equation  simplifies  for  propaga¬ 
tion  in  the  0  *  and  90  *  directions  when  one  also  considers  small  number  of  layers. 
In  this  paper  we  have  considered  a  five  layered  plate:  three  reinforced  layers  of  the 
same  thickness  and  two  interface  isotropic  soft  equal  thickness  layers.  Results  ue 
presented  using  both  an  analytical  technique  and  the  stiffness  method. 


Fig.  1.  Geometry  of  the  laminated  plate  with  interface  layers,  m  is  the  total  number 
of  layers  (m>5). 


2.  GOVERNING  EQUATIONS 


We  consider  time  harmonic  wave  propagation  along  the  x  or  y  axis.  Because  of  sym¬ 
metry  the  problem  reduces  to  two  imcoupled  ones:  plane  strain  in  which  the  dis¬ 
placement  components  are  u^,  o,  u„  and  SH  when  the  only  nonzero  displeurement 
component  is  u,.  In  this  paper  we  consider  only  the  plane  strain  case. 

2.1  Stiffness  Method 


Consider  the  ith  lamina  bounded  by  s=Z|_i  and  z=z,. 
this  lamina  will  be  given  by 
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The  stress  strain  relation  in 
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where  e-y  and  Cy  are  the  stress  and  strain  components,  respectively,  and  we  have  writ 
ten  7„=2e,a.  Note  that  if  the  y  axis  is  the  axis  of  symmetry,  then 


Then  the  problem  is  equivalent  to  that  of  an  isotropic  one.  In  order  to  get  good 
numerical  results  each  lamina  will  be  divided  into  several  sublayers,  p,  say.  Within 
the  ith  sublayer  we  will  choose  a  local  coordinate  with  the  origin  at  the  mid-plane 
and  Xj,  xj,  sj,  parallel  to  the  global  x,y,z  axes,  respectively.  Let  2hj  be  the  thickness  of 


this  sublayer.  Denoting  to  be  the  displacement  at  a  point  in  the  jth  lamina  we 
write 


Ci)  r  ^  f  ^  1  f  J.  1  cpj?  ^  , 


where  f„  (n=l,...,4)  are  cubic  polynomials  in  the  local  coordinate  zj  given  by 
fi  =  -^  (2-3fjj  +  tif).  fj  =  -^  (2+3»jj  -  fj-») 

fs  =  ■^  (l-»7j  -  »7j*  +  »?/),  ^4  =  •^  (“1  -^j  +  »?’  +  nf) 

Here  »?j=Zj/h  j  and  uj,  wj,  Xj.  Oj  are  the  values  of  Ujj,  u,,  a„,  and  at  the  jth  node. 
These  nodal  values  of  the  displacement  and  traction  components  are  functions  of  xj 
(=x)  and  t.  In  this  paper  it  will  be  assumed  that  the  time  dependence  is  of  the  form 
w  being  the  circular  frequency.  The  factor  will  be  dropped  in  the  sequel. 

The  equations  governing  the  nodal  generalized  coordinates  {uj,  Xjt  are 

obtained  using  Hamilton’s  principle.  The  details  can  be  found  in  [3].  It  is  found 
that,  if  {Q}  is  the  vector  of  all  the  generalized  coordinates,  and  if  we  assume  for  {Q} 

the  form 

{Q}={Q  >«"“  (3) 

then  {Q  }  satisfies  the  equation 

(k‘[K  I  -  ac>|Kj  -  k>[K  I  +  ik(K  I  +  (KJ)  {Q  }  =  0  (4) 

Matrices  [K  ],  etc.,  have  been  defined  in  (3).  [K  ],  (K  ],  and  [K  ]  are  linear  in  w*. 

For  nontrivial  solution  the  determinant  of  the  coefficient  matrbc  must  be  zero. 
The  solution  to  this  equation  provides  the  dispersion  relation  between  k  and  w. 


2.2  Exact  Solution 

Consider  wave  propagation  in  one  of  the  layers  along  one  of  the  principal  directions 
as  described  above.  Then  we  may  assume  that 

uM  =  (5) 

u  •'>  =  U, 

Further,  if  we  let 

=  A'e'*’*'  +  (6) 

Of  =  A'e""**'  -  A' e“‘*‘ 


(6) 


(7) 


=  B'e'“*  +  B'‘e“* 
nf  =  B'e'“^'  -  B''e“* 


then  we  may  write 

Ui  =  Vi'nr  +  vrnj- 

U3  =  v,'nf +  v,”nf 


Here  the  constants  Vj',  Vj”,  Vj',  and  Vj"  may  be  taken  as 

k|-k-a-s* 


Vj'  =  ik.  V3'  = 


isjj 


V,”  =  -is,.  V3“  =  - 


k|-k*a-3* 


ik5 


(8) 


S|  and  s|  are  given  by 

=  •il-(kVk.’(l+/S))±((kVk.’(l+«)’-‘W(k.’-'‘’)W-k>«))‘'"l  (9) 

In  the  above  we  have  defined 

a  =  4i/c^’  0  =  c^/cjS’  k,’=pwVcl? 

T=l+a0-^,  4=l+c|j//c  ^  . 

Note  that  if  the  y»axis  is  the  symmetry  axis  then 

=  cis/.  =  cl?  “  2cij) . 

In  that  case  the  solution  is  the  same  as  that  for  an  isotropic  medium. 

Using  the  solution  given  by  equations  (7)  and  (6)  and  applying  the  appropriate 
boundary  coniUtions  of  continuity  of  tractions  and  displacements  at  the  interfaces 
between  the  layers,  and  the  traction-free  boimdary  conditions  at  the  free  surfaces  of 
the  plate,  one  arrives  at  the  dispersion  equation  governing  k  and  u.  For  the  sake  of 
brevity  details  will  be  omitted  here. 


3.  NUMERICAL  RESULTS  AND  DISCUSSION 

As  an  application  of  the  techniques  described  above,  we  consider  a  fiber-reinforced 
plate  when  the  fibers  are  aligned  with  the  x-direction  (0*).  Properties  of  the  plate 
are  given  in  Table  1. 


Table  1,  Properties  of  0*  and  90*  laminae,  and  the  interface  layer.  All  the 
stiffnesses  are  in  the  units  of  10“N/m 


Cll 

cn 

Cl2 

C44 

Css 

0*  lamina 

1.2 

1.6073 

0.1392 

0.0644 

0.0350 

0.0707 

Interface 

1.8 

0.0865 

0.0865 

0.0475 

0.0195 

0.0195 

90*  lamina 

1.2 

0.1392 

1.6073 

0.0644 

0.0707 

0.0350 

Figure  2  shows  the  real  and  imaginary  branches  of  the  frequency  spectrum  for 
propagation  in  the  x-direction.  Solid  lines  in  this  figure  are  the  results  of  the  numeri¬ 
cal  method  and  the  open  circles  and  x’s  are  the  exact  solutions.  Note  that  the  two 
agree  very  well.  It  is  seen  that  as  k-*oo  the  slopes  of  the  first  symmetric  and 
antisymmetric  modes  tend  to  the  ratio  V^/  where  V  ^  is  the  Rayleigh  wave 

velocity  in  the  x-direction.  Results  by  the  stiffness  method  were  obtained  using  15 
sublayers. 


Fig.  2.  Dispersion  curves  for  the  real  and  imaginary  wave  numbers,  x:  symmetric 

kH 

and  o:  antisymmetric  modes.  n= - .  ■ '  ■■  '  ■  and  7=-;^,  H  being  the  thickness  of 

2>rv/(c5s//»Jo- 


the  plate. 


To  examine  the  effect  of  the  interface  soft  layers  we  then  consider  a  plate  com¬ 
posed  of  three  equal  thickness  unidirectional  (0  * )  fiber-reinforced  layers  separated  by 
two  thin  soft  layers.  Properties  of  each  layer  are  as  in  Table  1.  Figures  3a-3c  show 
the  dependence  of  the  phase  velocity  in  the  x-direction  on  the  frequency  when  there 
are  no  interface  layers  and  when  the  ratio  of  the  thicknesses  of  an  interface  layer  and 
a  lamina  takes  two  values,  0.1  and  0.2.  It  is  seen  that  the  presence  of  the  soft  inter¬ 
face  layers  b  to  lower  the  cut-off  frequencies  and  the  phase  velocities.  Abo  it  lowers 
the  surface  wave  velocity.  The  results  shown  in  Fig.  3  are  obtained  using  the  exact 
solution. 


Phase  Velocity 


Symmetric 
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Symmttnc 


Normalized  Frequency 

Fig.  3b.  Variation  of  c  with  -^kjH  when  the  interface  layer  thickness  ratio  is  0.2. 

It 

Using  the  stiffness  method  we  also  studied  the  dispersion  of  guided  waves  in  a 
0  /90  /O*  laminated  plate  with  interface  layers.  Fig.  4  shows  the  frequency-wave 
number  dependence  when  the  interface  layer  thickness  is  0.1.  These  curves  are  quite 
different  from  those  shown  in  Fig.  2.  To  examine  the  behavior  of  the  first  three 

branches  in  more  detail  we  show  in  Fig.  5  n=s — ■- — —  vs.  r*—-  It  is 

mterestmg  to  observe  the  slowing  down  of  the  waves  as  the  wavelength  becomes  of 
the  order  of  the  interface  layer  thickness. 


^7 
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Fig.  5.  First  three  real  branches  of  the  frequency  spectrum  for  a  0  *  /90  *  /O  *  plate 
with  soft  thin  interface  layers. 


4.  CONCLUSION 

Dispersion  of  guided  waves  in  a  fiber-reinforced  laminated  plate  with  or  without  soft 
interface  layers  between  the  laminae'^fiffl  beett'studied  «her^  using  an  exact  method  of 
solution  and  an  approximate  stiffness  method.  It  is  shown  that  the  approximate 
numerical  method  provides  a  solution  that  agrees  very  well  with  the  exact  solution. 
The  advantage  of  the  approximate  method  is  that  it  is  easily  used  for  any  number  of 
laminae  with  arbitrary  properties.  The  results  for  a  laminated  plate  with  soft  inter- 
f^e  layers  show  significant  lowering  of  cut-off  frequencies  and  phase  velocities.  Velo¬ 
city  of  surface^wav^is  found  to  decrease  with  increasmg  interface  layer  thickness. 
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